Influence of the particle number on the spin dynamics of ultracold atoms 
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We study the dependency of the quantum spin dynamics on the particle number in a system of ultracold 
spin-1 atoms within the single-spatial-mode approximation. We find, for all strengths of the spin-dependent 
interaction, convergence towards the mean-field dynamics in the thermodynamic limit. The convergence is, 
however, particularly slow when the spin-changing coUisional energy and the quadratic Zeeman energy are 
equal, i. e. deviations between quantum and mean-field spin dynamics may be extremely large under these 
conditions. Our estimates show, that quantum corrections to the mean-field dynamics may play a relevant role 
in experiments with spinor Bose-Einstein condensates. This is especially the case in the regime of few atoms, 
which may be accessible in optical lattices. Here, spin dynamics is modulated by a beat note at large magnetic 
fields due to the significant influence of correlated many-body spin states. 
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I. INTRODUCTION 

The spin degree of freedom of spinor Bose-Einstein con- 
densates (BECs) alows one to explore magnetism of ultra- 
cold quantum fluids [1], which exhibit intriguing phenom- 
ena, like the formation of coreless vortices [2] and other 
spin textures [3-5], and spontaneous symmetry breaking af- 
ter a quench [6]. Beside the investigation of the ground-state 
phases [7-10], the main focus of research is on the study of 
spin dynamics [11-15]. So far most experiments confirmed 
the mean-field (MF) description, including a nonlinear res- 
onance phenomenon near a critical magnetic field [16-19], 
which is caused by the interplay of spin-changing collisions 
and quadratic Zeeman shift. Interestingly, the Hamiltonian 
of the spin-dependent interatomic interactions appears also in 
nonlinear quantum optics [9]. As a result, spinor BECs allow 
one to explore, e. g., four-wave mixing [18, 20] and parametric 
down conversion [21-23] with matter waves. Moreover, spin- 
changing collisions are responsible for the dynamical evolu- 
tion of squeezed collective spin states and entanglement from 
uncorrected (product) states [24-28], which provides a way 
to overcome the standart quantum limit in precision measure- 
ments with matter waves. Those beyond-mean-field correla- 
tions can be described by means of the Bogoliubov approxi- 
mation as long as the quantum fluctuations of the spinor field 
operator are small [21-23, 29]. 

Correlations limit the validity of the MF approximation. 
Hence, knowledge of the boundaries of the Gross-Pitaevskii 
equation (GPE) is of great interest. The validity of the GPE 
has been proven for weakly interacting spinless bosons in the 
limit N —i' oo with Na fixed [30, 31], where N is the num- 
ber of particles and a is the scattering length. We show in 
this article by means of a numerically exact diagonalization of 
the effective spin Hamiltonian [9, 10, 29, 32], that the quan- 
tum spin dynamics in the single-mode approximation (SMA) 
converges towards the MF dynamics in the thermodynamic 



limit (TDL). This is, interestingly, the case for all strengths of 
the spin-dependent interaction (within the SMA). The conver- 
gence is, however, particularly slow in a regime, where the 
spin-changing collisional energy and the quadratic Zeeman 
energy are equal. We determine the validity time of the ini- 
tial MF dynamics under these conditions, which grows loga- 
rithmically with the number of particles. From this we expect 
that quantum corrections to the MF dynamics may play an im- 
portant role in experiments with spinor BECs. Additionally, 
we discuss a beat-note phenomenon in the spin dynamics of 
few atoms {N ^ 10), which may be observed in deep optical 
lattices. 

The paper is organized as follows. In Sec. II we outline the 
method to calculate the quantum spin dynamics. Afterwards, 
in Sec. Ill, we discuss few particles. First, in Subsec. Ill A we 
apply the formalism to two particles to make the method clear 
and to discuss similarities with the MF dynamics. Then, we 
discuss the spin dynamics of three particles in Subsec. IIIB, 
which is modulated by a beat note at large magnetic fields. 
A similar beat-note phenomenon is recovered for few atoms, 
which is shown in Subsec. IIIC. In Sec. IV we turn to the 
comparison of the A^-particle quantum dynamics with the MF 
dynamics. This is done for two typical initial states. First, in 
Subsec. IV A, we study the initial state, where all atoms are 
in the to = Zeeman sublevel, and then, in Subsec. IV B, we 
analyze the initial dynamics of the transversely magnetized 
state. We finally summarize our results in Sec. V. 



II. METHODS 

A. Effective spin Hamiltonian 

The two-body interaction between ultracold spin-1 atoms is 
modeled by a spin-dependent 6 potential [7] 

Vintin - r2) = 5{ri -~ f2)[hco + /1C2/1 • ^) 
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with the interaction strengths cq and C2 and the dimensionless 
spin-1 matrices fi of atom i ~ 1,2. Typically, cq is one or two 
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orders of magnitude larger than C2 . The atoms are confined by 
a spin-independent trapping potential Va-ipir)- Additionally, a 
homogeneous magnetic field along the z-direction generates 
the potential 

Vz = -hpf,~hq{A~f^), 

where p (x B and q oc sue the coefficients of the linear 
and quadratic Zeeman energy, respectively. The many-body 
Hamiltonian is then 



H" = 



Zm 



i 

Yhc26{n - rj)fi ■ fj- 



i<3 



i<j 



The first line contains the spin-independent part of H", which 
acts only in position space, and the second line contains the 
Zeeman Hamiltonian, which acts only in spin space. Only the 
spin-dependent interaction in the third line couples the spin to 
the motional degrees of freedom. 

In the absence of the spin-dependent interaction (c2 — 0) 
the states of the ground-state multiple! are of the form 

ipoiri, . . . ,rAr) (g) \xs) 

with ipo being the totally symmetric (nondegenerate) ground 
state of the spinless problem and with \xs) being an arbitrary 
totally symmetric iV-particle spin function [33]. In many ex- 
perimental situations one can restrict the description to these 
states [14, 15, 18, 19, 34, 35], since |c2| < cq. As a conse- 
quence, the motion of the atoms is frozen in the ground state 
and the system is essentially zero-dimensional. 
Let us assume that we have found the spatial ground state 
4>Q and the corresponding energy Eq . An integration over the 
spatial degrees of freedom leads to the effective spin Hamil- 
tonian (the diagonal offset Eq — ANhq is neglected) 



with 



C2 / drdfs . . .rN\'ipoir,r,r3, . ■ . ,rN)\ 
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The integral is the averaged local pair correlation function g'^^ 
of the ground state ipo- It can be viewed as the inverse volume 
of the ground state g^"^^ = 1/^- 

Using the projection operators Nm = X^i which 
count the number of particles with polarization m = +, 0, — , 
and the relation 



fN)"^ is the square of the total spin. 



where = (/i + . . 
we obtain a more convenient form of the Hamiltonian 



H' = -hp{N+ - N^) + hq{N+ + iV_) 
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{F^ - 27V). 



The z-component of the total spin F^ — — commutes 
with the occupation number operators Nm and the Hamilto- 
nian H', which leads to a decomposition of the dynamics into 
subspaces with different F^ — M — ~N, . . . ,N,i. e. 



{Nr, 



E(^" 



IM- 



M 



Therefore, the linear Zeeman energy, which is a constant in 
each subspace, does not influence the dynamics. In the fol- 
lowing we neglect the linear Zeeman energy and set p — Q. 
The remaining Hamiltonian 



H = hq{N+ + N_) + —{F^- 2N 



(1) 



has a spin-flip symmetry 

Hm 



M, 



since F^ is not changed by a rotation of 180° around the x- 
axis and since N+ + N- is unaffected if -r-> N- . 



B. Matrix representation 

The matrix elements of the Hamiltonian (1) are most con- 
veniently calculated within the second quantization formal- 
ism. We express the totally symmetric A^-particle spin func- 
tions \xs) by linear combinations of occupation number basis 
states |iV+, A^Oi which are eigenstates of the occupation 
number operators N„i — al^a^. The bosonic creation and 
annihilation operators aj„, act on these states in the usual 
way and obey the commutation relations [a^, aj^^,] = Smm' 
and zero else. 

Before we proceed, we change the labeling of the occupa- 
tion number basis states \N^, N^, N^) in order to simplify 
the following formulas. We use the set of quantum numbers 
(?7, M, N) instead of {N+, Nq,N^). Both labels are related 
to each other through tj ^ + M = - N_ and 
iV = + A^o + We will not explicitly refer to the 
number of particles A^, i. e. \N+, Nq, N^) = \rj, M). 

The first summand of (1), the quadratic Zeeman Hamilto- 
nian Hq — hq{N+ + A^_), is diagonal in the occupation num- 
ber basis and given by 

{r^,M\Hq\r^,M) ^T^hq. (2) 

The second summand, the spin-dependent interaction Hamil- 
tonian Hs = ^ - 27V' 
late we use the formula 



tonian Hs = ^ (F^ — 2 A^) , is tridiagonal. In order to calcu 



F^ = Ft 
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F+F^ +F_F+ 
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with the angular momentum creation and annihilation opera- 
tors F± = Fx ± iFy, which are given by 



F± = V2 (a±ao + ^o^q 



in dimensionless units. The diagonal elements of Hs are 

(77, M\Hs\v, M) ^ ^ [Af2 - 2if + f]{2N - 1)] (3) 

and the secondary diagonal elements are 

{7^ + 2,M\H,\t],M) = {t],M\H,\7^ + 2,M) = 

^y/{N -f]- 1){N - 7j){rj + AI + 2)(r? - M + 2) . (4) 

As discussed before, H decomposes into subblocks, which 
can be diagonalized independently for each eigenvalue of the 
total magnetization AI, since H commutes with F^. The basis 
states of one subblock are given by 

||M|,M),||M| + 2,A./),...,|r,„,„M) (5) 

with ?;max = — 1 or A^, leading to a total dimension of 

N - \AI\ 



dim(iV, AI) 



with the common floor function. 



1 



(6) 



C. Population dynamics 

Once we have determined the eigenstates | e) and eigenfre- 
quencies of H, we can calculate the time evolution of the 
system, if the initial state Itpi) has been specified. The initial 
state evolves according to 

\^Jj{t)) =ejip{^iHt/h)\^,), 

where the spectral representation of the time evolution opera- 
tor is given by 



exp(— iiJi/Zi) = yj |e)(e| exp(— iwji). 



(7) 



The time evolution of the relative population = Nq /N is 
thus given by 



(8) 



The population of the other spin components is completely de- 
termined by the conservation of the total number of particles 
N and of the total magnetization (F^) via 



n±{t) 



l±{F,)/N~no{t) 



For the initial states chosen here, (F^) = 0. In the following 
discussion we present only the time evolution of no and we 
neglect any contributions, which are constant in time, since 
they are determined by the initial state \tpi). 



D. Initial states 

We will discuss the time evolution of the number state 
I^Af) = |0, iV, 0), where all the atoms are in the m = Zee- 
man state, and the transversely magnetized state 



1 

7M 



2 + 



1 



N 



+ -al\ |0,0,0), 



where all the spins are pointing into the positive x-direction. 
This state is a superposition of number states from all sub- 
spaces 



^max 



N 

M=-N ri=\M\Mr)=2 



with the coefficients 




(9) 



E. Dimensionless coupling parameter 

The essential parameter, which characterizes the interplay 
between the quadratic Zeeman energy and the spin-dependent 
interaction energy, is given by 

K = > — ^ for large 

{2N-l)gs C2P ^ ' 

where p = N/V is the particle density. K can be positive or 
negative depending on the sign of q and C2 . 



III. FEW-ATOM DYNAMICS 
A. Two atoms 

We begin with two atoms to illustrate the method and 
since typical features of the two-atom dynamics occur also for 
larger particle numbers and in the MF limit. Two-atom spin 
dynamics was experimentally investigated in Refs. [14, 15] 
for a system being initially in the number state. 

Let us start with the calculation of the matrix elements of 
the Hamiltonian (1). The matrix H decomposes into 5 subma- 
trices with total magnetization AI — 2, 1, 0, —1, —2. Accord- 
ing to Eq. (6) the subspace with AI = has dimension 2 and 
the others have dimension 1. If one writes down Eq. (8) for 
an arbitrary initial state, one sees, that only energy differences 
within the same subspace lead to a sinusoidal oscillation with 
frequency ujij ~ — ujj\. Thus, only the AI = Q subspace 
contributes to the time evolution with exactly one frequency 
uj and the others lead to constant offset amplitudes, which we 
will neglect in the following. According to Eq. (5) the two ba- 
sis states of the M = subspace are |0,0) and |2,0). Using 
Eqs. (2-4) we obtain the 2x2 matrix 







^9s "^q-g, 



(10) 
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coupling parameter \K\ 

FIG. 1: (color online). Two-atom amplitude as a function of \K\ for 
the initially prepared number state, K < (red) and K > (blue). 




coupling parameter \K\ 



FIG. 2: (color online). Absolute value of the two-atom amplitude \A\ 
as a function of for the transversely magnetized state, K < 
(red) and K > (blue). Inset: comparison with the MF limit. 



Let US consider an arbitrary initial state, which is given by 



= a|0,0) +/3|2,0) + other terms. 



(11) 



(The other terms are the components of the irrelevant sub- 
spaces, like 7|1, 1) + d\2, 2) + . . . and so on. Thus, in general 
jap + ^ 1.) From the diagonahzation of the matrix (10) 
we obtain the spectral representation of the time evolution op- 
erator (7), which we insert into Eq. (8) together with the initial 
state (11). The result is an ordinary cosine oscillation 



nQ{t) — Acos{ujt) 



(12) 



with the amplitude 

2 



A= ^{2y25R(a/3*)(7<?,+.g2[2(| 



, ....... -\m-V2n{af3*)]l 

(13) 

where 5R(7) denotes the real part of 7 and 7* is its complex 
conjugate, and the frequency 



+ 1 



Time-independent terms have been neglected in Eq. (12), i. e. 
n'o(t) = no(t) - no,const.- 

Let us first discuss the initial state \62) = |0, 2, 0) = |0, 0), 
i. e. [a = 1, P — 0). In that case, the system undergoes Rabi 
oscillations [14, 15] and the amplitude (13) becomes 



converges to zero on the one hand, since the number state is 
an eigenstate of Hq, and on the other hand, since the occupa- 
tion number operator no(t) = a^a^^/N commutes with Hq, 
and thus becomes a constant of motion. 

For the transversely magnetized state IC2) we calculate the 
coefficients (a = 1/2,(3 = •\/2/4) using Eq. (9). The ampli- 
tude (13) takes the form 



A 



qgs 



{2q 



K 



A is negative for K < 0, which corresponds to a phase shift of 
TT in the cosine function of Eq. (12), and is not in contradiction 
to the requirement no(t) > 0, since tt-o, const. > 1^1- 

Fig. 2 shows 1^1 as a function of |if | for K < (red) 
and A' > (blue). The transversely magnetized state is an 
eigenstate of the spin-dependent interaction Hg and thus the 
oscillation amplitude is small for small \K\. Again, the ampli- 
tude drops down in the opposite limit of large \K\, since the 
(]^4) occupation number operator no{t) — aQa^/N commutes with 
Hq and becomes a constant of motion. Between these limit- 



ing cases the amplitude has a maximum, which is located at 
I A' I = 1. The MF amplitude (App. A) has a maximum at the 
same position and shows the same limiting behavior, which 
is clear, since the above arguments were independent of the 
number of particles (inset of Fig. 2). 



A 



(2g-.g,)2 + 8.g2 9 _ 6if + 9A'2 ' 



(15) 



Fig. 1 shows the amplitude A as a function of the coupling 
parameter \K\ for the two cases AT < (red) and AT > 
(blue). For small \K\, the Hamiltonian is dominated by the 
spin-dependent interaction Hg = ^ (F^ _ 2A^) , and since 
the number state is not an eigenstate of Hg, the system under- 
goes large oscillations. In the opposite limit of large \K\, the 
Hamiltonian is approximately equal to the quadratic Zeeman 
energy Hq — hq{N^ + iV_). Here, the oscillation amplitude 



B. Three atoms 

In the case of three atoms, there are 3 subspaces of dimen- 
sion 2 (those with M — 0,±1), which contribute to the dy- 
namics. However, since H is symmetric under polarization 
reflection M o —M, the M = ±1 subspaces yield the same 
frequency. The corresponding Hamiltonians are given by 



i/0 = fi 



f V6gs 
\V6gs 2q + g 



, H^ 







2gs 2q - 3gs 



2gs 
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FIG. 3: (color online). Evolution of three atoms in the transversely 
magnetized state. (a)ylo/i as a function of A', (b) ljo/i asafunction 
of K. (c) Dynamics for = —5 (red) compared to the MF evolution 
(blue). The oscillation of the m = Q population is modulated by a 
beat frequency of 2gs, which is absent in the MF limit (App. A). 



where in we have subtracted the offset h{q + 2gs) from 
the diagonal. For the initially prepared number state, the dy- 
namics is restricted to the M ~ {) subspace, which leads to 
similar results as in the previously discussed two-atom case. 

By contrast, the transversely magnetized state jCs) is dis- 
tributed over all subspaces. Since all subspaces evolve inde- 
pendently, they can be evaluated in the same way as before 
and added thereafter. The time evolution of the m = popu- 
lation is now given by a sum of two cosines 

7iQ(t) — Aq cos(woi) + Ai cos(cjii) 

(constant offsets are neglected) with different amplitudes 

and unequal frequencies 

c^o = v/(2g + .9.)2 + 24g2, = ^(2? - 3.g,)2 + 16.g2 . 

The amplitudes Aq and Ai as a function of K show the same 
behavior as discussed for two atoms [see Fig. 3(a)]. The fre- 
quencies and ijji as a function of K are shown in Fig. 3(b). 
At large K they differ by 4gs so that the oscillation dynamics 
is modulated by a beat frequency of 2173: 



^{cos[(2g + g,)t] + cos[(2g - ig,)t] } 



1 
5^ 



cos[(2g - gs)t\ cos{2gst) . 



(16) 




0.2 0.3 

time t [27r/|.g.,|] 

FIG. 4: (color online). Evolution of seven atoms in the transversely 
magnetized state for — —5. The exact numerical calculation 
(blue) is compared to the first-order perturbative evolution (17) (red). 
The envelope function of (17) is drawn as a dashed line. 



C. Dynamics of few atoms for large K 

For intermediate values of K, the dynamics of few atoms 
looks rather chaotic, since many frequencies contribute to the 
evolution. For large K, the transversely magnetized initial 
state shows a fast oscillation, which is modulated by a beat 
frequency, as in the three-atom case. The limiting dynamics 
is obtained from a perturbative calculation: For 1 /K « 0, the 
quadratic Zeeman energy is the dominant part of the Hamilto- 
nian and thus we choose the number states |iV+, A^o, N^) as a 
basis. The interaction energy is a small perturbation. In a first 
step we approximate the eigenstates and energies up to first 
order in A = N/{8K). This is separately done for subspaces 
of H with different magnetization = M, since their evo- 
lution is decoupled. The result and the representation of the 
initial state (9) is inserted into Eq. (8). After neglecting 
terms of order , A'^ , . . . one obtains the time evolution 



iV- 1 



2K{2N - 1) 



cos[(2g - gs)t\ [cos{2gst)\ 



N-2 



Fig. 3(c) shows the three-atom dynamics (red) compared to 
the MF dynamics (blue) at K = —5. 



(17) 

The amplitude of the oscillation is proportional to 1 /K, as 
in the two-particle and MF limiting cases (Fig. 2). The fre- 
quency of the fast oscillation, which is approximately given 
by « 2g, is determined by the level spacing of Hg, which co- 
incides again with the two-particle (14) and MF results (A3). 
But different from these two limiting cases, the fast oscillation 
is modulated by a beat frequency of 2gs. 

Fig. 4 shows the time evolution of seven atoms at K ~ —5. 
The beat note is a clear signature that coiTelated spin states 
contribute to the dynamics of few atoms. We believe, that the 
beat frequency can be observed in deep optical lattices with 
few atoms at each lattice site, similar to the measurements of 
Refs. [14, 15]. The observation is facilitated by the fact, that 
the beat frequency is independent of the number of particles. 

Eq. (17) is valid for rather large particle numbers. We found 
excellent agreement with the initial evolution of < 300 parti- 
cles for K — ~h from a comparison with the numerical re- 
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FIG. 5: (color online), (top) Evolution of the number state with TV = 
10, 50, 200 bosons at K = 0. The time scale is 2n/\gs\. (bottom) 
Evolution for TV = 250, 500, 1000, 2000 bosons at isT = 0. The time 
scale is now Ttdl = 27r/(|3s|TV) = 2tv/{\c2\p). In the MF limit, 
the number state is a steady state, and no(t) = 1. With increasing 
particle number TV the population no(t) decreases more slowly on 
the time scale Ttdl until it stays constant in the thermodynamic limit. 



suits, although A 1 for these parameters. 



IV. TV-PARTICLE QUANTUM DYNAMICS VS. 
MEAN-FIELD DYNAMICS 



One can show by means of the large-TV method [36] that the 
MF description of the spin system considered here becomes 
exact in the thermodynamic limit (TDL) [37]. That means 
that the time evolution of a coherent initial state is described 
by the MF equations of motion in the TDL. Further, one of the 
ground states of the system becomes a coherent state and its 
energy can be calculated by minimizing the expectation value 
of H on the MF phase space. Product states (Al) become 
identical with coherent states in the limit TV — > oo and thus the 
TV-particle quantum dynamics of the initial states considered 
here converges towards the MF solution in the TDL. In the 
following we study the influence of the finite particle number 
TV on the quantum corrections to the MF dynamics. 



A. Number state 

The number state is a steady state of the MF equations of 
motion for all values of the coupling strength K (App. A). 
That makes it particularly useful to study corrections, which 
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FIG. 6: (color online). Relative population of m — ±1 Zeeman 
states, n± = 1 — no, after t = 0.5 Ttdl as a function of {~K). 
Beyond-MF corrections are largest around K = —1 and decrease 
with increasing particle number. 



go beyond the conventional MF dynamics, since quantum spin 
fluctuations are strongly amplified in this initial state [22, 23, 
29]. 



1. Time evolution for zero K 

At K = 0, the Hamiltonian consists only of the interac- 
tion Hs, which has the total-spin states \F, M) as eigenstates. 
To calculate the TV-particle quantum dynamics, one needs the 
representation of the occupation number operator ?io and the 
initial state \6m) in the total-spin basis {\F, T\f)}. The calcu- 
lation is done in App. B. The evolution is given by 



no 



{t) = Y,CFC0s[g,{2F + i)t], 



where the coefficients Cf are approximated by 

„ F f IF^' 
Op 



— cxp , 
TV ^ V 2 TV 



(18) 



(19) 



for large particle numbers TV (see App. B). The dynamics is 
periodic with 2ii/\gs\, since all the frequencies in (18) are 
multiples of Qs. Fig. 5{top) shows a rapid oscillation on the 
time scale 2'n/\gs\, which seems to be in contradiction with 
the MF result, since the number state is a steady state of the 
MF equations of motion for all K. The time scale 27r/|gs|, 
however, becomes infinite in the TDL, since 2'kV/\c2\ oo 
for V ^ oo. Thus, we need a time scale, which stays constant 
in the TDL. A proper choice is given by 

TTDL = 2V(l5s|iV)-2^/(|c2|p), 

which is inversely proportional to the spin-dependent inter- 
action energy. By plotting the initial TV-particle quantum dy- 
namics as a function of Ttdl [Fig- 5(bottom)] one sees conver- 
gence towards the MF limiting behavior with increasing parti- 
cle number TV. The numerical result is confirmed by the ana- 
lytical formulas: The distribution (19) has a maximum around 
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FIG. 7: (color online). Evolution of the m = population for 2000 
particles in the transversely magnetized state for different coupling 
strengths: K = -0.1 (red), K = -1 (green) and K = -2 (blue). 
The corresponding MF evolution is drawn as a dashed line. 



F = VN and a width proportional to v N. Thus, the only 
cosine oscillation, which contributes to the dynamics (18) has 
zero frequency in the TDL, since 

2VNgs = 2c2p/Vn (for N oo). 



2. Time evolution for large K 

A perturbative calculation with the small parameter A = 
1 / {2K) leads to the first-order result 



4(Ar- 1) 
K^{2N -I) 



2qt 



2Af- 3 
K{2N - 1) 



As in the case of two particles the amplitude drops down pro- 
portional to \/K^ [compare with Eq. (15)] and the frequency 
of the oscillation converges to 2q. Moreover, the amplitude is 
proportional to and thus becomes zero in the TDL. 



3. Beyond-mean-field corrections in dependence of K 

We finally analyze the connections to the MF dynamics in 
dependence of the coupling parameter K. For that reason, 
we have numerically calculated the fraction of atoms in the 
TO = ±1 Zeeman states, n± (t) = 1 — rio(i), after a fixed time 
t — O.STtdl as a function of {—K). Fig. 6 shows a resonant 
enhancement of the to = ±1 population around K ^ —1. For 
N = 1000 the relative population of the to = ±1 Zeeman 
states is largest and the resonance maximum is slightly below 
\K\ = 1. With increasing particle number N the resonance 
maximum decreases and its position converges to K = —1. 

Beyond-MF coiTections were analyzed in a similar way in 
a recent experiment [22, 23] and explained within the Bo- 
goliubov theory. Different from here, two resonances have 
been observed there, which is due to the restricted motion of 
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FIG. 8: (color online). Oscillation period of the initial evolution in 
dependence of {—K). For small and large \K\, the A^-particle and 
MF period are in good agreement. In the vicinity of K — —1, the 
period is strongly enhanced in systems with large particle numbers. 
The enhancement of the oscillation period is absent in small systems 
with less than 20 particles. 



the atoms to the ground state in our approach. However, on 
the first resonance, the coiTesponding Bogoliubov mode has 
a similar shape than the MF ground state so that our method 
is applicable there. Indeed, in a large trap, both approaches 
lead to the same value for the position of the first resonance, 
namely |g/(c2p)| = 1. 



B. Transversely magnetized state 

1. Time evolution for three different values of K 

The transversely magnetized state shows a rich dynamics 
in the MF limit [18, 34]. Fig. 7 shows the initial evolution of 
N = 2000 atoms for three different coupling strengths. We 
compare the exact iV-particle quantum dynamics (solid line) 
to the corresponding MF evolution (dashed line). 

At small \K\ (red) one sees no difference between the two 
solutions within the time 5 Tjdl- The evolution is an ordinary 
cosine oscillation with amplitude A = K/A and frequency 
u! = 2c2P period Tjdi^/2), which is determined by the 
spin-dependent interaction energy (App. A). 

In the opposite limit of large \K\, the MF evolution (blue 
dashed) is given by a cosine oscillation with amplitude A ~ 
1/(47^) and frequency w = 2g period Ttoi^/{2K)]. The 
exact A'^-particle quantum dynamics (blue solid) oscillates 
with the same fast frequency of 2q. Moreover, the fast os- 
cillation is modulated by the envelope function of Eq. (17). 
The beat note vanishes in the TDL, since 



cos 



N-2 



1 



2{c2ptf 
N 



Despite the beat note, both oscillations coincide quite well for 
times, which are smaller than 2 Ttdl- 
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FIG. 9: (color online). Initial evolution of the relative population 
of the m = Zeeman state on resonance (K = — 1) for different 
particle numbers A'^ compared to the MF limit. 
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FIG. 10: (color online). Position of the first oscillation maximum 
against the number of particles A'' on resonance | = 1). 



At K = —1, the relative population of the m — Zeeman 
state converges towards 1 in the MF limit (green dashed). The 
amplitude becomes maximal with A — 1/4 and the oscillation 
period diverges. The period of the exact TV-particle quantum 
evolution (green soUd) is enhanced, but still finite. Moreover, 
the green solid curve does not converge towards 1, but rather 
oscillates around 0.82. Note, that the evolution is coherent and 
exhibits a revival although the initial dynamics seems to be 
strongly damped and to converge towards a constant value. At 
K — —1, one observes the largest deviations from the MF dy- 
namics. The initial evolutions coincide only for times, which 
are smaller than Ttdl/2. A similar observation was made in 
the experiment [18, 19]. 



paiticle numbers smaller than « 20, the maximum is ab- 
sent. One fuither sees that the position of the resonance max- 
imum rapidly approaches 1 from below with increasing parti- 
cle number TV. 

The enhancement of the oscillation period has been exper- 
imentally observed in an antiferromagnetic spin-2 ^^Rb BEC 
[18] and, for a similar initial state, in a spin-1 ^''Na BEC [19]. 
The large differences between the exact -particle quantum 
dynamics and the MF evolution close to the critical coupling 
strength at |iir| = 1 due to finite- coiTections may be rele- 
vant for studies of quantum chaos [38, 39]. 

3. Limit for the validity of the mean-field dynamics 



2. Beyond-mean-field corrections in dependence of K 

Fig. 8 shows plots of the oscillation period against the cou- 
pling strength {—K) for different paiticle numbers and the 
MF limit. In the MF limit, the oscillation period was obtained 
from the analytical solution (App. A), while in the other cases, 
the period is 2 times the position of the first maximum of the 
initial oscillation. 

One sees that the A^-particle oscillation period coincides 
quite well with the MF period in the limit of small and large 
\K\, for all particle numbers. In the resonance region around 
K = —1, the deviations are quite large for small particle num- 
bers; see the curves for A^ = 2, 20 and 50. With increasing 
particle number, the A/-particle period coincides almost ev- 
erywhere with the MF limit, except in a small region around 
K = —1, which shrinks to zero in the TDL; see the curves 
for A^ = 300 and 2000. However, exactly on resonance at 
K = —1, the deviation from the MF limit is always infinitely 
large for finite particle numbers A^. 

The enhancement of the oscillation period is a typical fea- 
ture of the spin systems with large particle numbers. One 
clearly sees the first occurence of a weak maximum for A^ = 
20, which becomes rather pronounced for N — 2000. For 



We conclude from the previous discussion, that the conver- 
gence to the MF limiting dynamics in systems with a finite 
number of particles is slowest at |_fir| = 1; see Figs. 6-8. 
Hence, we study now, how fast the oscillation period con- 
verges towards infinity at the critical value of \K\ — 1, since 
this worst-case scenario provides us with the smallest upper 
bound for the validity of the MF approximation. 

Fig. 9 shows the initial evolution of the relative population 
of the m Zeeman state at = -1 for A^ = 50 - 2000 
paiticles. As expected, the initial evolution coincides with the 
MF dynamics for longer times, when the paiticle number is 
increased. One sees that the A^-paiticle quantum dynamics is 
close to the MF evolution for times, which are smaller than 
the position of the first oscillation maximum. 

Hence, we plot in Fig. 10 the position of the first oscilla- 
tion maximum against the number of particles on a double 
logarithmic scale. One finds a logarithmic dependency of the 
validity time Tk{N) of the MF approximation on the number 
of particles 

TK=-i(iV) = Ttdl [0.13 + 0.09 logioW] 
T;^=i(iV) = Ttdl [0.18 + 0.09 logioW]. 

Thus, for a spin-1 BEC of 10^ atoms and a negative coupling 
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parameter, K < 0, one finds the MF dynamics to be valid for 
times, which are smaller than w 0.67 Ttdl- 



V. CONCLUSIONS AND OUTLOOK 

We studied the quantum spin dynamics of ultracold spin-1 
atoms within the single-mode approximation. This was done 
on the basis of an exact diagonalization of the effective many- 
body spin Hamiltonian [9]. The chosen method allowed us 
to discuss the crossover from few atoms to small condensates 
within one framework. 

Our numerical calculations showed convergence of the 
quantum spin dynamics towards the mean-field dynamics 
in the thermodynamic limit for all strengths of the spin- 
dependent interaction. Moreover we showed that quantum 
corrections to the mean-field dynamics are particularly large 
at the critical value of the coupling parameter, where the spin- 
changing collisional energy equals the quadratic Zeeman en- 
ergy. For the state, where all the atoms are in the m = sub- 
level, we compared the initial quantum spin dynamics to re- 
sults obtained from a Bogoliubov approximation [22, 23, 29]. 
For the transversely magnetized state we estimated the valid- 
ity time of the initial mean-field dynamics at the critical value 
of the coupling parameter, which grows logarithmically with 
the number of particles. From this estimate we conclude that 
quantum corrections to the mean-field dynamics may play an 
important role in experiments with spinor BECs. It would be 
interesting to study systematically the dependency of the spin 
dynamics on the particle number in future experiments. 

Quantum corrections to the mean-field spin dynamics are 
particularly large in the regime of few atoms. Here, many- 
body spin correlations lead to a beat-note phenomenon at large 
magnetic fields. The regime of few atoms {N ^ 10) should 
be accessible in deep optical lattices. 

Our results may be relevant for the study of quantum chaos. 
It was shown [38] that the strong nonlinear behavior of the 
mean-field equations at the critical value of the coupling pa- 
rameter leads to classical chaos in the spin dynamics of spin-2 
atoms. But exactly at this point, the finite-iV quantum dynam- 
ics largely deviates from the classical limit. It would be inter- 
esting to study this aspect systematically, similar to the work 
in Ref. [39] on a periodically driven double-well system. 
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Appendix A: Mean-field dynamics 

One assumes in the MF approach, that the system is in a 
product state of the form 



(Al) 



with a,n being complex numbers, which are normalized ac- 



cording to J2 r, 



1. For the states (Al) and the Hamil- 



tonian (1), one derives the MF equations of motion [34] 

idta^ = gsN(^A*aQ + (f^) a_^_^ + qa^ 
idta^ = gsN(^Aaj^ + A*a_^ (A2) 
idta_ = gsN(^Aa^- {f^) + qa_ , 

where we have defined A = (/+)/\/2 = {a*^aQ + agCt-) 
and (/,)^(|a+|2-|a_p). 

The number state |0,iV, 0) = jO)**^ is a steady state of 
(A2), since 

ao(i) = l, a±(i)=0 

is a solution of the MF equations of motion (A2) for all values 
of gs and q. Thus, this state shows no population dynamics in 
the MF limit and na{t) = |ao(i)P = 1- 

The time evolution of the transversely magnetized state 



|Q = (-l/2|+)-i/^/2|0) + l/2|-)) 

is calculated in Refs. [34, 35]. In this state, all the spins are 
pointing into the positive y-direction. The solution of the MF 
Eqs. (A2) is given in terms of Jacobi elliptic functions [40] 



ao(i) = 



V2 



cnfc(f)dnfc(f) i(l + fc)snfc(f) 
(l-fc)snfc(f) 



l + fcsn2(f) 
icn,(f)dn,(f; 



l + fcsn^(f 



where s ~ exp(— i(gsiV — q)t/2) and k = 1/K. For the spin 
populations n„i{t) = \am{t)\'^ the solution simplifies to 



no{t) = [l-ksnUqt)]/2, 
n±(t) = [l + ksnl{qt)]/4. 



(A3) 



These solutions are also valid for the initial state |Cjv), which 
is used here, since the evolution of the relative populations 
is unaffected by rotations around the z-axis. 

For small k ^ 1/K, one can approximate snfc(a;) « sin(x) 
and Eq. (A3) becomes 



no{t) 



1 

Ik 



1 

Ik 



cos{2qt) (large if). 



That means, the evolution is a cosine oscillation with ampli- 
tude A — 1/ {4:K) and frequency uj = 2q. 
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For large k = 1/K, we approximate snfc(a;) « sm{kx)/k, 
which leads to 



no{t) 



1 K 

2 " T 



K 



cos(2c2pi) (small if). 



That means, in the interaction dominated regime, the oscilla- 
tion amplitude is A = K/A and the frequency is uj = 2c2p. 

At \K\ = 1, the evolution becomes aperiodic and the rel- 
ative m — population asymptotically approaches 1, i. e. 
no(t) — 1 for t cx). Here, the dynamics exhibits a maxi- 
mum of the amplitude and the oscillation period diverges. 



Appendix B: Quantum dynamics of the number state at zero K 

In the following, we derive the population dynamics of the 
number state at zero K. In this limiting regime, the Hamilto- 
nian consists only of the interaction H^, which has the eigen- 
basis {\F, A/)}. The occupation number operator A'o and the 
initial state \9n) need to be expressed in this basis to calculate 
the dynamics. The expansion of \9n) is given by 



N 



\eM) = |o,iv,o) 



E xf\f,q). 



(Bl) 



F=F„,i„,A_F=2 



Since \9m) has the eigenvalue M ~ 0, it is a superposition 
of the states \F,M = 0). Due to symmetry reasons, the sum- 
mation runs over F = Fmin, i^,ni„ + 2, . . . , with F^i^ = 
or 1 if A^ is even or odd, respectively [33]. The coefficients 
Xf are determined later. By inserting the expansion (Bl) into 
Eq. (8) the evolution becomes 

Mt) = ExFXF'(i^,0|A^o|i^',0)cos[(u;i.-wj^Oi]/^ 

F,F' 

(B2) 

with the frequencies ojp — gs\F{F + 1) — 2A^] . The matrix 
elements {F, 0|A"o|-F''7 0) are calculated in the following. 

The occupation number operator Aq can be written in terms 
of spherical tensor operators 



From the Wigner-Eckart theorem one finds, that the matrix 
elements {F' ,M'\tI^^\f,M) vanish for |F - F'\ > 2. The 
nonzero matrix elements of are 



(F,0|r(j^^|F,0) and {F + 2,0\T^''>\F,0) . 
The Wigner-Eckart theorem allows one to calculate the matrix 

(2) 

elements of Tq ' from a special class of matrix elements 



(2) 



{F + q,0\4^^\F,0) 



{F,2;0,0\F + q,0) 



{F,2;^F,-q\F + q,~F-q) 



x{F + q,-F-q\T[^^\F,-F) 



(B5) 



for q = 0,2. The brackets {fi, f2',''ni,m2\f' ,m') are the 
Clebsch-Gordan coefficients (CGC). Eq. (B5) simplifies the 
calculation, since the states \F, —F) have a rather simple rep- 
resentation in the occupation number basis 



\F,-F)=CF{alf[alal-{alY/2 
with the normaUzation constant 



|0,0,0) (B6) 



E 

fe=0 



JV-F 



■fc!(fc + F)!(A^-i^-2fc)! 

2JV-F-2fe ■ 



(B7) 



After inserting Eqs. (B3), (B4), (B6) and the required CGCs 
into Eq. (B5), a lengthy calculation leads to the nonzero ma- 
trix elements 



W0|Tr|f.O> = -L(3f-2iV + 6|)^ 



(B8) 



and 



{F + 2,0\4^^\F,0) 



3 CF N~F (F + l)(F + 2) 
2cp+2 2 y {2F + l){2F + 3)' 



(B9) 



(2) 

where Tq is the zeroth component of the one-particle spher- 

(2) 

ical tensor operator Tq ' of rank 2. Its five components are 



t!j^2 = '^Va 



^±1 - y^iflo'^T ^ "±"0 



1 



7^(2) _ 



76 



(B3) 



(B4) 



A^ is proportional to the identity matrix and thus its matrix 
elements are 



(i^,0|A^|F',0) = N5_ 



FF' 



where 



fc=0 



N-F \ 2 
2 

k 



' k\{k + Fy.{N - F - 2k)\ 

2N-F-2k 



(BIO) 



The coefficients xf of the expansion (Bl) can be calculated 
by means of Eq. (B6). Using 



\F,0)^l/^/{2Fy.{F+y\F,-F) 



(BID 



we obtain 



XF = (i^,0|0,7V,0) = \0,N,0)/./{2Fy. 



CF 



1 



2F 
F 



(B12) 
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FIG. 11: (color online). Exact amplitudes Cf [(B13), blue crosses] 
of 2000 atoms compared to the approximation [(B24), red line]. 



The evolution is obtained by inserting Eqs. (B8)-(B12) into 
Eq. (B2). Only terms with \F — F'\ =2 lead to non-constant 
contributions to the evolution: 



JV-2 



no{t)' 



CFCOs[gs{2F + 3)t], 



F=F„,i„,A_F=2 

with frequencies gs{'2F + 3) = — i^f and amplitudes 



(7V-F)4(7V-1)!(-) 



xJryT^'y'/^^T^'\- (bi3) 

\\Fj \F + 2j {2F + l){2F + 3) 

The amplitudes (B13) can be approximated for large N. 
The most involved part is the approximation of Ci?. The fac- 
torials in (B7) can be written as binomial coefficients. 



kl{k + F)l{N -F- 2k)\ ^ N\ {j^^^ (^,% 
where i9 = 2fc + F. We obtain from (B7) 



N 

E 



iV! 



F i}=F,\-d=2 ^ 2 



N-F\ 2 
2 



2 



(B14) 

The terms with d N dominate the sum in (B14). The ap- 
proximation 

1 {k-n/2f' 
2 



1 (n 
2" U 



1 



27rv/ri/4 



exp 



n/4 



(B15) 



holds for large n. Applying this to {(^^p-^j^^ gives 



1 ^ 



TV! (F'^ 



(B16) 

Further, we approximate d = N m the exponential and the 
square root of (B16). The exponential dependence on F^ 



shows, that c^? is negligible for F N. Thus, we assume 
F <C in the following. The binomials in (B16) are ex- 



panded into factorials and approximated using {n 
nl/n'' for n k: 



- fc)! 



^-F 
2 



N 



2N 



We apply (B15) to the right-hand side of (B17) 



Further, for F N, one gets to first order 

'N-F\^-^ ( F{N-^) 



\ N ) 



cxp 



N 



(B17) 



(B18) 



(B19) 



We insert the approximations (B17)-(B19) into (B16) and ap- 
proximate the sum by an integral 



.^F F^/{2N)^'^ 

4 ^ 2~ 2 



N 



N 



e " , 



where an additional factor 1/2 accounts for the step size of 2 
in the sum. The integral has the value 



N 



N 



FJL 
6 ^ 



F{N - 1?) 



(B20) 

with the error function erf(a::). Inserting the upper limit N of 
the integral into (B20) leads to erf(O) = 0. With F N, the 
lower limit F gives erf{VF) w 1 for F > 5 and thus 



7f ^"pUaT 



2^ 2 



(B21) 



The maximum of (B21) is at F = ^J~Nj2. The Gaussian 
exp(— F-^/A^) has a width of ^ N/2, which is an upper bound 
for the width of cp- Now we approximate the other terms in 
(B13). Since F « ^/N for the dominant contributions. 



2F\ /2F + 4 
F / I F + 2 



2F 
F 



i2F 



(B22) 



where we used (B15) in the last step. For 1 <C F <C A^, the 
remaining factors become 



A^-F / (F+l)(F + 2) 
N ' Y (2F+ l)(2F + 3) 

Inserting (B21HB23) into (B13) leads to 

F 



N 



exp 



'2 A^ 



(B23) 



(B24) 



The approximation (B24) and the exact amplitude (B13) are 
compared for 2000 atoms in Fig. 11. 
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